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We study a recently proposed measure for the quantification of quantum non-Markovianity in the 
dynamics of open systems which is based on the exchange of information between the open system 
and its environment. This measure relates the degree of memory effects to certain optimal initial 
state pairs featuring a maximal flow of information from the environment back to the open system. 
We rigorously prove that the states of these optimal pairs must lie on the boundary of the space of 
physical states and that they must be orthogonal. This implies that quantum memory effects are 
maximal for states which are initially distinguishable with certainty, having a maximal information 
content. Moreover, we construct an explicit example which demonstrates that optimal quantum 
states need not be pure states. 
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I. INTRODUCTION 

Many quantum physical systems are modelled approx- 
imately as closed systems. However, in most practical 
cases the interaction of a given system with its envi- 
ronment can not be neglected and this makes it neces- 
sary to use the theory of open quantum systems [ ] . In 
the past, it has been very common to neglect quantum 
memory effects and to resort to a Markovian approxi- 
mation. In many applications the time evolution of the 
open system is described by a quantum dynamical semi- 
group represented by a generator of the Lindblad-Gorini- 
Kossakowski-Sudarshan form [2, 3]. 

However, when the environment has a nontrivial struc- 
ture, one has to account for memory effects influenc- 
ing the dynamics of the open system. During recent 
years, there has been significant conceptual, theoretical, 
and experimental progress dealing with non-Markovian 
processes [4-20]. In particular, the very definition of 
non-Markovianity and quantification of quantum mem- 
ory effects in the dynamics of open systems has received 
a lot of interest. Several measures for quantum non- 
Markovianity have been proposed which are based on 
different mathematical and physical concepts [6, 8-11]. 

In the present paper we investigate the non- 
Markovianity measure proposed in [ ] which expresses 
the degree of memory effects in terms of the amount of 
information exchanged between the open system and its 
environment. Within this approach the degree of non- 
Markovianity is connected to certain optimal initial pairs 
of quantum states, leading to a maximal flow of informa- 
tion from the environment to the open system. Here, we 
will study the mathematical and physical properties of 
such optimal state pairs. 
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Throughout this paper % refers to the Hilbert space 
of the open quantum system and the corresponding set 
of physical states, i.e. the set of positive trace class oper- 
ators p with unit trace is denoted by S(H). We assume 
that the dynamics of the open quantum system can be 
represented in terms of a one-parameter family 

* = {$t | < t < T} (1) 

of completely positive and trace preserving (CPT) linear 
maps $ t , where <J>o is equal to the identity map [ ]. We 
will use the quantity 

A/"($)= max / dta(t, Pl ,p 2 ) (2) 
Pi,a€S(n)J a>0 

as a measure for the degree of memory effects in the open 
system dynamics [8, 9]. Here, 

Cr(t, Pl , P2 ) = -V(<f>M,$ t ( P2 )), (3) 

where T> denotes the trace distance [21, 22] defined by 

V( Pl ,p 2 ) = ~Tr\ Pl - P2 \. (4) 

By definition, the non-Markovianity measure A/'( ( f ) ) is a 
positive functional of the process $. It is zero if and only 
if the trace distance between any pair of initial states is 
a monotonically decreasing function of time t signifying 
a continuous loss of information from the open system to 
its environment. A nonzero value for the measure means 
that there is an initial state pair for which the trace dis- 
tance increases over a certain time interval which can be 
interpreted as a flow of information from the environ- 
ment back to the open system implying the presence of 
memory effects. The time integral in (2) determines the 
total backflow of information for a certain pair of initial 
states p\ and p 2 . The quantity J\f(Q) is then found by 
taking the maximum over all initial state pairs, i.e. A/"($) 
represents the maximal possible backflow of information. 
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A pair of states p%,p2 € S(W) is said to be an opti- 
mal state pair if the maximum in Eq. (2) is attained for 
this pair. In the present paper we are interested in the 
mathematical and physical features of the maximization 
procedure involved in expression (2). In particular, we 
will obtain general statements about the physical proper- 
ties of optimal state pairs which give rise to the maximal 
possible degree of memory effects in their dynamics. To 
this end, we first show in Sec. II A that optimal pairs of 
states must lie on the boundary of the state space S(H). 
We then proceed in Sec. II B to demonstrate that the 
states of any optimal pair must be orthogonal. This is 
our central result which is physically very plausible since 
it implies that an optimal pair must have the largest pos- 
sible initial trace distance T>(pi, p%) = 1. Thus, the max- 
imal flow of information from the environment back to 
the open system emerges if the initial state pair is dis- 
tinguishable with certainty, i.e. has a maximal informa- 
tion content. We discuss specific features of the case of 
infinite-dimensional state spaces in Sec. II C, and present 
in Sec. II D an alternative proof for the orthogonality of 
optimal state pairs which employs the joint translatabil- 
ity of nonorthogonal states. 

In the simple case of a qubit, the orthogonality of the 
optimal pair implies that both states of the pair must be 
pure. In Sec. Ill we study the question of whether this 
statement holds true for higher dimensional systems. The 
answer is negative. In fact, we will construct in this sec- 
tion an example for the dynamics of a three-level system 
(A-system) for which the optimal pair is not a pure state 
pair. Thus, optimal state pairs exhibiting maximal back- 
flow of information can indeed be mixed states. Finally, 
we draw some conclusions from our results in Sec. IV. 



II. CHARACTERIZING OPTIMAL PAIRS OF 
QUANTUM STATES 

A. Restriction to the boundary of the state space 

The state space S(T-l) of an open quantum system is 
given by the set of positive trace class operators on its 
Hilbert space W with trace one [23], 

S(H) = {p | p > 0, Trp = 1} . (5) 

On the basis of the convex structure of the state space 
one can define the boundary dS(H) of S(H) as follows. 
A point p £ S(H) is defined to be an interior point if and 
only if for all a £ S(H) there is a real number A > 1 such 
that 

(1- X)a + Xp€S(H). (6) 

Denoting the set of all interior points by S(H), one de- 
fines the boundary of the state space by 

dS(H) = S{H)\S{H). (7) 



Thus, p £ dS(H) if and only if there exists a £ S(H) 
such that for all A > 1 the operator (1 — X)cr + Xp does 
not belong to the state space S(^H). 

It is easy to show that a density matrix p 6 S(H) 
belongs to the boundary dS(W) of the state space if it 
has a zero eigenvalue. In fact, let \ip) be a normalized 
eigenvector of p with eigenvalue zero, and P = \tp){ip\ the 
corresponding projection operator. Then, for all A > 1 
the operator (1 — X)P + Xp has the negative eigenvalue 
1 — A < with corresponding eigenvector \ip), and thus 
does not belong to the state space. 

For finite-dimensional systems the boundary dS(W) 
actually contains only states with zero eigenvalues. This 
can be seen by setting A = 1 + e with e > in (6) which 
yields the operator p' = (1 + e)(p — e(l + e) _1 cr). If 
all eigenvalues of p are strictly positive, it follows from 
the continuity of the roots of the characteristic poly- 
nomial that for sufficiently small e also all eigenvalues 
of p' are positive. This shows that states with strictly 
positive eigenvalues are in the interior S(H) such that 
dS(H) = {p E S(H)\0 £ spec(p)} for finite-dimensional 
systems. 

As has been mentioned already in the introduction a 
pair of states pi , pi is said to be an optimal state pair if 
the maximum in Eq. (2) is attained for this pair, i.e., if 

= [ dt a(t,p 1 ,p 2 ) (8) 

holds for this pair. If we speak of an optimal pair we will 
always assume that the quantum process under consid- 
eration is non-Mar kovian, i.e. Af($) > 0, since otherwise 
any pair of quantum states would be trivially optimal. 
Moreover, we note that the optimality property implies 
of course that p\ ^ P2- 

Theorem 1. Let pi,P2 G be an optimal state pair. 

Then both states lie on the boundary of the states space, 
p 1 ,p 2 £dS('H). 

Proof. Suppose that (at least) one state of the pair, say 
P2, does not belong to the boundary. Hence, p2 is an 
interior point and there exists A > 1 such that 

p 3 = (1 - X)pi + Xp 2 € S{H) (9) 

is a quantum state (see Fig. 1). The time evolution of 
the three states is given by pi{t) = § t (pi), i = 1,2, 3. By 
the linearity of the dynamical map we have 

P3 (t) = (1 - X) Pl (t) + X P2 (t) (10) 

and, hence, 

p 1 (t)-p 3 (t) = X(p 1 (t)-p 2 (t)). (11) 

It follows that 

V(p 1 (t),p 3 (t)) = XV( Pl (t),p2(t)). (12) 

Note that A is a fixed number strictly larger than 1. Thus, 
the last equation tells us that the trace distance between 
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pi(t) and Pd,{t) is always larger by the constant factor 
A than the trace distance between pi(t) and p 2 (t). This 
implies that the quantity J a>0 dt o~(t, pi, p^) is larger than 
J a>0 dt a(t,pi,p2) by the same factor A (see Eq. (3)). It 
follows that p\ , pi cannot be an optimal pair which is a 
contradiction. Consequently, any optimal pair of states 
must belong to the boundary of the state space. □ 




Figure 1. Illustration of the decomposition (9) and of the 
corresponding time evolution given by Eq. (10). 

Thus we see that the maximization over all initial state 
pairs in the definition (2) of the non-Markovianity mea- 
sure Af(<&) can be restricted to the boundary dS(T-l) of 
the state space, 

7V($) = max / dt a(t,p 1 ,p 2 ). (13) 

We emphasize that this result is valid for any Hilbcrt 
space T-L and for any family of linear dynamical maps $ t . 
The linearity guarantees the invariance of the decompo- 
sition defined in (9) for all times t (see Fig. 1). 

B. Orthogonality of optimal state pairs 

In this section we demonstrate that optimal state 
pairs must be orthogonal which strengthens the result 
of Sec. II A. We recall that two quantum states p\ and 
P2 are said to be orthogonal, p\ _L p 2 , if and only if their 
supports are orthogonal, supp(pi) _L supp(p2), where the 
support is defined as the subspace spanned by the eigen- 
vectors with nonzero eigenvalues. Note that the orthog- 
onality of two states implies that both states have a zero 
eigenvalue and, hence, that both states belong to the 
boundary of the state space. We also mention that the 
trace distance between any two quantum states satisfies 
T>(pi,p 2 ) < 1, where the equality sign holds if and only 
if p\ _L p2 ■ Orthogonality is thus equivalent to unit trace 
distance [21, 22]. 

Theorem 2. Optimal state pairs are orthogonal. 

Proof. Let p\,p 2 € S{T-L) be an optimal pair of states, 
and suppose that the states are not orthogonal, pi JL p 2 . 
According to the Jordan-Hahn decomposition [22] there 



exist positive and orthogonal operators Pi and P 2 such 
that 

Pi-p2=Pi-P 2 . (14) 

Taking the trace of this equation we see that A = TrPi = 
TrP2- From p\ ^ p 2 we conclude that A > 0. Using 
Eq. (14), the orthogonality of Pi and P 2 , and the fact 
that the trace distance between nonorthogonal states is 
always strictly smaller than 1, we find 

1 > V{ Pl ,p 2 ) = V{P U P 2 ) = X - (TrPr+TrPa) = A. 

(15) 

Thus, we have < A < 1. Now we define the operators 
(Ji = Pi /A and a 2 — P 2 /X. Being positive and of unit 
trace, these operators represent quantum states. More- 
over we have 

CTl - 0"2 = j{pi - p 2 )- (16) 

By use of the linearity of the dynamical maps we obtain 

a 1 (t)-a 2 (t) = j(p 1 (t)-p 2 (t)), (17) 

from which it follows that 

V{a 1 (t),o- 2 (t)) = ~V(p 1 (t),p 2 (t)). (18) 

Since A -1 > 1 we can conclude from the last equation 
that the pair a\ , a 2 yields a non-Markovianity which is 
strictly larger than that of the pair p\ , p 2 , which con- 
tradicts the assumption that pi,p 2 is an optimal pair. 
Hence, pi and p 2 must be orthogonal. □ 

It follows from theorem 2 that the maximization in 
the definition (2) of the measure Af($>) for quantum non- 
Markovianity can be restricted to orthogonal initial state 
pairs. Again, this result holds for any Hilbert space H 
and any family of linear dynamical maps $ t . 

C. Infinite dimensional state spaces 

It is important to emphasize that for infinite dimen- 
sional spaces an optimal state pair need not necessarily 
exist. To be mathematically more precise, the maximum 
in the definition (2) for the non-Markovianity measure 
should be replaced by the supremum to account for such 
cases. 

In theorems 1 and 2 we have assumed that an op- 
timal pair of states does exist. If there is no such 
pair there is a sequence of state pairs Pi,p 2 such that 
J a>0 dt a(t, pi, p 2 ) converges to A/"($) for n — > oo. Em- 
ploying the construction used in the proof of theorem 2 
one can show that the pairs Pi,p 2 can always be taken 
to be orthogonal. This means that the non-Markovianity 
measure can be approximated with arbitrary precision by 
orthogonal state pairs and that we can write 

JV($)= sup j dt a(t, Pl ,p 2 ). (19) 

P1-I-P2 Jcr>Q 
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D. Orthogonality and parallel translations 

In this section we present an alternative proof for the 
orthogonality of optimal state pairs, where we restrict 
ourselves to finite dimensional Hilbert spaces, dimH = 
N. This proof could be of interest also in other contexts 
since it relies on the behavior of pairs of states under 
parallel translations. 

The idea of the proof is based on the observation gained 
from two-level systems that nonorthogonal states on the 
boundary of the state space can be simultaneously trans- 
lated by a traceless Hermitian operator to yield a pair of 
mixed states, while the trace distance of the pair is in- 
variant under such translations. To make this idea more 
precise we define 

£(H) = {A | A^O, A = A\ TrA = 0} (20) 

to be the set of nonzero, Hermitian and traceless opera- 
tors on H. 

Definition. Two states p\,p2 € S(H) are called jointly 
translatable if and only if there exists an A € £{T~L) such 
that p k - Ae S(H) for k = 1,2. 

Hence, two states are said to be jointly translatable if 
and only if there is a nontrivial Hermitian and traceless 
operator A which can be subtracted from the states with- 
out leaving the state space. We prove that any pair of 
nonorthogonal states is jointly translatable in such a way 
that both translated states do not belong to the bound- 
ary of the state space. 

Theorem 3. If pi,P2 € S{T-i) and p\ JL p 2 then p\,p 2 
are jointly translatable. Moreover, there exists an opera- 
tor A e £{%) such that p k - A (£ dS{H) for k = 1, 2. 

Proof. Let pi j2 be given in terms of their spectral decom- 
position, 



»=i 



(21) 



where p\ , i = 1,...,N, denotes the eigenvalues and 

\ipj, k>> ) the corresponding eigenvectors of pk- The assump- 
tion p\ JL p2 implies that p\ and p% have at least one 
eigenvector with nonzero eigenvalue which are not or- 
thogonal. After possible relabeling one can assume that 

a = (ipi~^ 7^ and p^\p^ ^ holds. Moreover, 
by an appropriate choice of the phases of the eigenstates 
one can assume without restriction that a is real and 
positive, i.e. < a < 1. We consider the following super- 
positions of the two overlapping eigenvectors, 



iv>±> - c ± (i^) ± i4 2) >) 



(22) 



normalization constants obey c ± = 
a). The projections onto these states are de- 



where the 

noted by P±. By positivity of the overlap a one has 



c + < c_. Now, we define 
A f = e-B, 



P_ - 



N 



(23) 



where e > is a real number to be chosen later. Clearly, 
= B, B ^ and TiB = so that A e is a candidate 
for the traceless hermitian operator which simultaneously 
shifts both states. Furthermore, as c+ < c_ one recog- 
nizes that the two last terms are negative semidefinite 
while the first one is positive semidefinite. In order to 
show positivity of pk = Pk — A e for an appropriate choice 
of e we consider the quantity (x|/3fc|x) f° r an arbitrary 
normalized vector \\) € H. We find 



(XIMX) =B k (\x))+p[ k) \(x\4 k) ) 
+ e(^±)\ X \P_\ X ) + 



<x\p+\x) 

- m: 



N 

p[ k) \(x\4 k) )\ 2 

Re((^ 1) | X ) ( X |^ 2) ))} 

- pf\(x\4 k) )\ 2 



= B k (\x)) 



+ 4cie 



{%-\(x\4 k) )\} 



(24) 



where B k (\x}) = Ym=2pT' I {xWt') |" > and since 
c^ 2 — cl 2 = 4a by definition of c±. In the last step one 
uses the fact that Re(z) < \z\ for all complex numbers z, 

and that KV'i'^lx)! < 1 f° r all normalized vectors by the 
Cauchy-Schwarz inequality. Consider now the quadratic 
function 



(fc)| 



(k) 2 , , 2 f a 
Pi x +4c + e \^ 



i 



(25) 



One can show that this function is strictly positive for 
k = 1, 2 if e satisfies 



< e < 



a 



(fc) 

mm p\ . 



(26) 



Nc% k=iX 

Thus, if we choose any e satisfying (26) we obtain 
g € (\ (xl^i^) I) > for all normalized vectors |%) e % 
and k — 1,2 which leads to (xl/3i,2|x) > 0- This demon- 
strates that pi t 2 is positive and has no zero eigenvalue, i.e. 
Pi, 2 € S{%) due to the characterization of the boundary 
in terms of the eigenvalues given in Sec. II A. In partic- 
ular, pi and p 2 are jointly translatable. □ 

We remark that the converse of the statement of theo- 
rem 3 holds, too. Thus, two states are jointly translatable 
if and only if they are not orthogonal. 

Theorem 2 can now be proven as follows. Let p\,p2 
be an optimal pair, and suppose that pi / p 2 . Then by 
theorem 3 there exists an operator A € £(H) such that 
Pk = Pk — A $l dS{%s)- Applying theorem 1 we conclude 
that the states pk are not optimal. Since p\—p2 = pi —f>2 
it follows that also the states pk are not optimal, which 
represents a contradiction. Therefore, the optimal pair 
has to be orthogonal. 
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III. PURITY OF OPTIMAL PAIRS 

As a simple application of theorem 2 one obtains the 
result that for all non-Markovian quantum processes of a 
two-dimensional system (qubit) the maximal backflow of 
information occurs for a pair of pure, orthogonal initial 
states, corresponding to antipodal points on the surface 
of the Bloch sphere. This follows immediately from the 
fact that for qubits the set of pure states is identical 
to the boundary dS(H) of the state space. For higher- 
dimensional systems this is no longer true, i.e., the set 
of pure states represents a proper subset of the bound- 
ary in this case. In this section we will construct an 
explicit example for an open system dynamics in a three- 
dimensional Hilbert space for which the optimal pair is 
not a pair of pure states. 

We consider a A-system which interacts with an off- 
resonant cavity field. The weak-coupling master equation 
of this model is given by [9] 



dt 



pit) = -iAi(t)[|o><a|,p(t)] - i\2(t)[\a}(a\,p(t)} 



+ Mt)[\b)(a\p(t)\a)(b\--{p(t),\a)(a\} 

+ Tii(t)[|c)<o|p(t)|a)<c|-i{p(t),|o)<o|}] ) (27) 

where \a) refers to the excited and |&),|c) to the two 
ground states. The coefficients an d 71.2 (t) are 

determined by the spectral density of the cavity field. 
Introducing the functions 



/(*) 



= -(r)i(t)+r>a(t))/2 -i^W+iaW) 



gi{t )= ( d sll (s)e-( Dl{s)+D2{s) y, 



(28) 
(29) 



where 



r t ,t 
D t (t) = I ds 7i {s), Li(t) = / dsXi{s), (30) 
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we find that the solution of the master equation yields 
the dynamical map 

(\f(t)\ 2 Paa f(t) P ab f(t) Pa 

&t(j>)=\ fWplb 9l(t)paa + Pbb Pbc 



\ f(t)*P* a 



Pbc 



92{t)Paa + Pc 



(31) 

The functions /, g\ and (72 have to obey the following 
relations which guarantee that $^ is trace preserving and 
completely positive, 



9i(t)+92(t) + \f(t)\ 2 = l, 
9xAt) > 0. 



(32) 
(33) 



We consider now, for simplicity, the case where the 
Lamb-shifts Xi are equal to zero, while the decay rates 7$ 
in the dissipator of the master equation (27) are chosen 
to obey one period of small oscillation, 

A1.2W = 0, 7i, 2 (t) = 0.03 ■ sin(i). (34) 



It is easy to check that for this choice the conditions (32) 
and (33) are satisfied. We have carried out numerical 
simulations, drawing random pairs of pure, orthogonal 
initial states and determining the corresponding increase 
of the trace distance for each initial pair. The results are 
shown in Fig. 2. We see from the figure that there is a 
finite gap between the maximal possible increase of the 
trace distance for pure, orthogonal initial pairs, and the 
increase of the trace distance corresponding to the initial 
pair 

Pi = \a)H P 2 = l(\b)(b\ + \c)(c\), (35) 

which consists of the excited state and the uniform mix- 
ture of the two ground states (indicated by the arrow in 
Fig. 2). Thus, this example clearly demonstrates that 
for Hilbert space dimensions larger than two, the opti- 
mal initial state pair can indeed contain a mixed state 
for certain non-Markovian dynamical maps [25]. 




0.02 0.04 



0.06 0.08 
total increase 



0.10 0.12 



Figure 2. (color online) Probability of the total increase of the 
trace distance for randomly drawn pairs of pure orthogonal 
states for one period of oscillation of the decay rates. The 
arrow indicates the value of the maximal increase of the trace 
distance for the pair of states given in Eq. (35) including a 
mixed state. The sample size of pairs of pure orthogonal states 
is equal to 10 . 



IV. CONCLUSION 

In this paper we have shown that optimal pairs of ini- 
tial states for non-Markovian quantum dynamics must be 
orthogonal which implies that these states can be distin- 
guished with probability 1 by a single measurement. Op- 
timal state pairs thus have the maximal possible amount 
of initial information and are therefore capable of emit- 
ting and reabsorbing the maximal amount of information 
during the non-Markovian dynamics. 

We emphasize that the proof of these statements only 
relies on the convexity of the state space and on the lin- 
earity of the dynamical maps. Our results are thus very 
general and can be applied to any quantum process de- 
scribable by a family of linear dynamical maps $ t in any 
Hilbert space. Strictly speaking, even the complete posi- 
tivity of the dynamics is not needed. In fact, it suffices to 
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assume that the maps <!> t are positive since even trace- 
preserving positive maps are contractions for the trace 
distance [24]. 

We have further demonstrated that for Hilbert spaces 
with dimensions of at least 3 optimal state pairs need not 
consist of pure states, in contrast to the case of a qubit 
where optimal pairs are always antipodal points on the 
Bloch sphere, and as such pure. The example constructed 
here leads to an optimal pair consisting of a pure and 
a mixed quantum state. We conjecture that in Hilbert 
spaces of dimension 4 and higher one can also construct 
quantum processes for which both states of the optimal 
pair are true mixtures. Finally, we mention that on the 
basis of the present results further, more specific state- 
ments could be proven if one assumes additional prop- 
erties of the quantum process, such as invariance under 



certain symmetry groups, or the existence of invariant 
states. In addition, the notion of joint translat ability pro- 
vides new insights into the structure of the state space 
and might serve for further applications. 
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